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Abstract 

Starting from Buchmiiller’s observation that a chromoelectric flux tube me¬ 
son will exhibit only the Thomas type spin-orbit interaction, we show that a 
model built upon the related assumption that a quark feels only a constant 
radial chromoelectric field in its rest frame implies a complete relativistic ef¬ 
fective Hamiltonian that can be written explicitly in terms of quark canonical 
variables. The model yields linear Regge trajectories and exhibits some sim¬ 
ilarities to scalar confinement, but with the advantage of being more closely 
linked to QCD. 


Typeset using REVTj^X 


1 


I. INTRODUCTION 


The nature of quark confinement is an old but still active problem. The transition 
to a linearly rising potential at large quark-antiquark separation is evident from lattice 
simulations. The associated “flux tube” held distribution has long been an attractive starting 
point to develop classical and semi-classical meson models [jl| and non-relativistic hybrid and 
glueball models as well as relativistic quantized models of meson dynamics . Buchmiiller 
@] has proposed that a hux tube which is purely chromoelectric in its rest (co-moving) frame 
must yield a spin-orbit interaction that is only of the kinematic (Thomas) type. The result 
follows classically, as we review in the next section, from the given fact that there is no 
magnetic held in the co-moving frame. 

In the present paper we show that this simple picture is far more powerful than originally 
thought. We show that when one asks what helds in the stationary laboratory frame are 
equivalent to those of a co-moving chromoelectric hux tube, one is led to a complete relativis¬ 
tic ehective theory. Section || contains a relativistic version of Buchmuller’s argument based 
on the Thomas equation. The complete ehective Hamiltonian in the Bethe-Salpeter form 
is given in Section 13. and the reduction and consequent relativistic corrections are worked 
out in Section |I^. In this section we also compare our result with the relativistic corrections 
expected from scalar and time component vector conhnement as well as those expected from 
QCD. The Regge structure of our model is discussed in Section 0. Our conclusions are 
summarized in Section B 


II. BUCHMULLER’S PICTURE OF THE SPIN-ORBIT INTERACTION IN QCD 

The classical covariant treatment of a spinning charge moving in electric and magnetic 
helds was hrst introduced by Thomas in 1927. This formalism automatically includes the 
“Thomas precession” ehect. The interaction energy of the spin magnetic moment is ^ 

Hso = -s ■ G, (2.1) 

m 

where m is mass of the particle and 

G = (|-l + -)B-(f-^)vxE. (2.2) 

Here, and in what follows, the particle charge has been absorbed into the helds {e.g. eB ^ 
B), and for simplicity we assume v-B = 0, which always will be true in the case we consider. 


2 




The gyromagnetic ratio g is expected to take the value g = 2, but we consider the general 
case. 

The above result can be re-expressed in terms of the quark co-moving fields B' and E' 
by the usual Lorentz transformation 


B|| = Bj|, (2.3) 

E|| = Ej|, (2.4) 

Bx = 7(B1 + vxE'), (2.5) 

Ex = 7 (E'^ - V X B'). (2.6) 


Substituting the transformed fields into ( | 2 . 2 | ), using v ■ B = 0 and 7 ^ = 1 — we obtain 

9 1 


G = — B' + - V X E'. 

27 7 -h 1 


(2.7) 


Buchmiiller’s observation ^ follows directly from (E3)- If the flux tube is pure chromo¬ 
electric in its rest (co-moving) frame, then B' = 0 , and the spin-orbit term energy is by (| 2 . 1 | ) 
and (|2.7|) 


Hso = 


m(7 -h 1) 


s ■ (v X E') 


( 2 , 8 ) 


For low velocity quarks we have 7 ~ 1, L ~ mr x v. With the assumption E' = —ar', 
where a is the magnitude of the electric field, (sometimes called the string tension in flux 
tube models), the spin orbit energy becomes 


H ~ 

^-‘SO — 


2rn?r 


s L. 


(2,9) 


This pure “Thomas” spin-orbit term arises in a natural way in the Buchmiiller picture. In 
QCD confinement the moving quark does not “cut” electric field lines as it moves since the 
electric field is carried along in a roughly straight flux tube connecting the quarks. 

For completeness we also mention the spin-orbit terms arising from other assumptions 
often made about the confinement potential. With a Lorentz scalar confinement there is no 
magnetic field in any frame, so only the purely kinematic Thomas term is present and the 
spin-orbit term ( p.9|) again results. For the time component vector confinement B = 0 in the 
lab frame, and E' = —af ~ E. The Lorentz transformations ( |2.5|) then yield B' = —yv x Eh 
By (p.l|) and (p.7|) we then find for a low particle velocity 




vector 


2m?r 


(^-l)s-L. 


( 2 . 10 ) 
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For g = 2, the above equation becomes 


H, 


vector 


2m?r 


s ■ L, 


( 2 , 11 ) 


which is of the same magnitude but opposite sign from the pure Thomas term. 

Discussing QCD in an electrodynamic context is not unreasonable in this model of quark 
dynamics given the assumptions we make about the physical helds. The particular flux tube 
held conhguration present depends crucially upon the nonabelian nature of QCD while the 
interaction of the helds with the quarks does not in this case. This is so because we assume 
that the helds agree with that of a hux tube, wherein only one component of the QCD held 
strength namely the radial chromoelectric held, is nonzero in a co-moving frame. Thus 
we may write the SU{3) held strength as in a 

frame that is moving with the same velocity as the quark. As the held tensor has only one 
non-zero component in this frame, we hnd that the Lorentz components of the QCD held 
strength trivially all commute with one another. 


iGa'g'(x'), = inGl,p,G%,T<^ oc = 0. (2.12) 

Because the held strengths commute in one frame, they must commute in all frames. 

[Gag{x'), G^,{x')] = G^x')] = 0. (2.13) 


Thus we may perform a gauge transformation upon G °'to rotate it into an abelian 
subalgebra of SU{3). That is, there exists a gauge transformation D(a;) such that the held 
strength can be rotated to lie in, say, the “three” direction of SU (3), D(a;)G“^,y(a;)T“D“^(a;) = 
G^yi x)T^, or that the held strength is purely in a U{1) subalgebra of SU{3). Choosing the 
held strength to he in the “three” direction of SU{3) partially hxes the gauge, leaving a 
residual U{1) gauge freedom. We then make the identihcation of the “electromagnetic” held 
strength as E^i, = G^,^. 


III. THE EFFECTIVE HAMILTONIAN 

We begin this section by establishing a relativistically valid four-vector potential in the 
laboratory frame that yields B' = 0 in the instantaneous rest frame of the quark. While the 
four-vector potential is not a gange invariant object, the fnll wave equation for the moving 
quark is gauge invariant provided that we perform a gange transformation on the quark wave 
function as well as on the fonr-vector potential. In the instantaneous rest frame the E' and 
B' helds at the qnark are assnmed to be those given in the Bnchmiiller pictnre. 
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B' = 0, 

E' = -at'. 


(3.1) 

(3.2) 


It is important to note that while the quark is acted upon by the chromodynamic helds, it is 
a source of the chromodynamic helds as well. The vector potential depends on the position 
and motion of the quark as well as upon the held point r 

A^ = Af,{r,t,rQ,pQ). (3.3) 

The curl of the vector potential A with respect to the held point is the magnetic held 


B(r,l,rQ,pQ) = V, X A(r,f,ro,po). 


(3.4) 


It is not the general magnetic held (|3.4|) that appears in a wave equation, but rather the 
operator expression 


B.j(re,(,pe) = 'i(pQ X A(ro,(,ro,pQ) + A(re,(,re,po) x pg) = V,, x A(rg, (, rg, pg), 

(3.5) 


which we assume to be the magnetic held acting on the quark, 

B(rg, t, rg, Pg) = Vg X E(rg). 


(3.6) 


We hnd the simplest such potential yielding this magnetic held by applying a Lorentz 
transformation to the co-moving potentials^ 


x4g = ar\ 
A! = 0, 


(3.7) 

(3.8) 


and using the relation between momentum and quark velocity. 


Pq = iQmwQ + A, (3.9) 

to obtain the potentials in the lab frame. The only quantities that enter into the dynamical 
equations are functions of the quark position and momentum so we may drop the subscripts 


^It is straightforward to include an additive constant in the expression for A'q which we neglect 
for simplicity. 
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Q on quantities that refer to the quark from now on as there is no more possibility of 
confusion. 


Ao(r', p) = aqr' = ar'^J[l + (p/(m + ar'))^], (3.10) 

(IT^ 

A(r , p) = aqvr' = —-—^p. (3.11) 

m + ar 

We assume that E' oc r' and E oc r so we must use the Lorentz transformed radius in the 
lab frame 


r' = — = rJl — (v ■ f)2. 
T-ll 


(3.12) 


To lowest relativistic order, this correction is only important in Aq. 

Buchmiiller’s argument predicts the Thomas spin-orbit long range interaction, but it is 
not a complete dynamical theory. In particular, an effective Hamiltonian is needed to solve 
for energy levels, transition rates, etc. To proceed in this direction we promote the lab frame 
four-vector potential 24^(r) (he. the four-vector potential at the quark position r) given in 
( p.lOl ) and (I3.11D , to a suitably symmetrized operator, and use it with the Salpeter equation! 
in the heavy-light limit P, 


A+[q; ■ (p — A(r)) -1- Pm + Ao(r)ll — EljA+T = 0. 


(3.13) 


As before, we have absorbed the charge into the potential. In the above equation (which is 
also known as the “no-pair” equation 0) we have 


Eo±Ho 

Hq = cx ■ p + Pm, 
Eq = Jp^ -h m^. 


(3.14) 

(3.15) 

(3.16) 


The angular momentum is given by 

J = r X p + ^cr. 


(3.17) 


^Use of the Salpeter equation avoids possible inconsistencies of the Dirac equation due to the 
Klein paradox. For purposes of deriving a non-relativistic reduction either equation is acceptable 


because both equations have the same nonrelativistic reduction |]T|. 
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where the canonical momentum is the sum of the mechanical quark momentum and the field 


momentum A, as given in eq. (|3.9| ). 
If we dehne 


A = AqI — q: ■ a, 


(3.18) 


where 11 is the unit matrix, we can put ( p.l3|) in normal Salpeter form (with E denoting the 
light quark energy) 


{Hq + A+AA+ - = 0. 

Since the wave function already has the positive energy projected, it follows that 


A_^ = A_ 


/(r) 

^(r) 


= 0 , 


and hence 


1^0 + m 


By using ( p.21| ) in (|3.19|) we obtain the two component form of Salpeter equation 

Hf{r) = A/(r), 

where 


H — Eq + (A_|_AA_|_)ii + (A_|_AA_|_)i2 


o- • P 

Eq + tu' 


(3.19) 


(3.20) 


(3.21) 


(3.22) 


(3.23) 


After a few lines of algebra, using (|3.14|) - (p.l6|) and (|3.18|) , we obtain the effective Hamilto¬ 
nian in the form 


H = Eq + Hp 


Hm, 


(3.24) 


where 


He = 
Hm = 


2 Eq l 
1 


2A, 


0 


{Eq m)Ao + (cr • p)Ao(cr ■ p) 

1 

{Eq + m){(T ■ A)(cr ■ p) 


Eo + m\' 


Eo + m 


+ (cr • p)(cr ■ A) 


(3.25) 

(3.26) 


No gauge condition has been imposed upon the (17(1)) vector potential so the effective 
Hamiltonian given in (|3.24|) - (|3.26|) is valid for any gauge. The eigenstates of this Hamiltonian 
in two different gauges will be related to each other by (17(1)) gauge transformations in the 
usual way. 
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IV. RELATIVISTIC CORRECTIONS 


To make contact with previous work we consider the case where the fermion has a large 
mass m and moves slowly. We wish to hnd an effective action that can be used perturbatively 
to compute relativistic corrections to the non-relativistic reduction as an expansion in the 
small quantities p/m and ar/m. 

The expansion in p/m and ar/m is facilitated by a few commutator identities. Given 
that 

rlT^ r) 

[p^F(r)] = -V2F(r)-2—(4.1) 
and since the expectation values satisfy 

2{~)=-{V^F), (4.2) 

dr or 

we hnd 

([p^F(r)]) = 0. (4.3) 

Hence, to lowest order in p^/m? ~ m^^(l ■Ep^/2m? ...)), 

(|i5„±‘,F(r)|)~0. (4.4) 


It is this relation that guarantees that the Dirac equation and the Salpeter equation have 
the same non-relativistic reduction. 

With these aids we can immediately reduce the terms of the effective Hamiltonian (|3.24|) , 




p 


p 


2m 8m^ ’ 
~ Eo + Ho ■ 


-(p-A +A-p) 


1 dAo 1 _ . 


(4.6) 

(4.6) 


8 m^ 2m' 

The potentials we found are functions of quark position and momenta, which is a situation 
not often encountered in electrodynamics. To construct a quantum theory with Hermitian 
potentials, we symmetrize the products of position and momentum. Ideally we would use 
Weyl ordering to make connection with the midpoint prescription for evaluating a path 
integral. Using ( p.lO| ) and ( |3.11| ), we hnd 


V X A(r) ~ — , 
mr 

Ho(r) ar — -— (r ■ p -|- p ■ r) 
4m 


^{-rir.PiPj+PiPjrir,-^ 


(4.7) 

(4.8) 










Am? 
a 


{j?r + rp^ 


A ~ -—(pr + rp) , 

2m 


V^An(r) ~ — . 


2a 

r 


(4.9) 

(4.10) 

(4.11) 


The second and third terms in Aq result from the 1/7|| correction to the radius. 


With these expressions the reduction ( |4.6[) becomes 


hf ~ Eq + 

~ Eo + a?" — 


3a 


arp^ 


arp 


2m^r 2m^ 2m^ 2m^r 

3a aL^ a 

-s ■ L. 


s • L, 


(4.12) 

(4.13) 


2 m^r 2m‘^r 2m‘^r 

For contrast, we compare these results to those obtained in scalar and time-component 
vector confinement and in QCD. 

If the conhnement is pure time component vector we substitute 


into (|4.6|) to obtain 




vector 


Aq = ar, A = 0, 




4m^r 2 m^r 
In the case of scalar confinement we take 


s L. 


A similar reduction gives 

hfgcalar — Eq -|- QT 


A = (3ar. 


arp^ 


aL? 


s L. 


(4.14) 


(4.15) 


(4.16) 


(4.17) 


4m^r 4m^ 4m^r 2m^r 
Finally, the low quark velocity expansion of the Wilson minimal area law provides the 
effective Hamiltonian implied by QCD ||11|| , 

a aL^ 


hfqCD ^ Eo + ar + 


s L. 


(4.18) 


36m‘^r Qm‘^r 2m^r 
We observe that our model, scalar confinement, and low quark velocity expansion of the 


Wilson loop model all have the pure Thomas spin-orbit interaction. The Darwin terms 
differ but are sensitive to the ordering prescription applied to symmetrize the operators. 
We note that even in the case where numerical coefficients do not agree, our model has 
coefficients of the same sign as both the scalar conhnement and Wilson loop results. 
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V. REGGE BEHAVIOR 


The slope of the leading Regge trajectory can be calculated classically in the limit of 
highly excited circular orbits as 

, L 


a = 


^2’ 


(5.1) 


where E is the excitation energy and L is the orbital angular momentum. For such orbits 
spin and quark mass are negligible and, by (|3.10 ) and ( 3.11 ), we have in that limit 

Aq —> + p2, (5.2) 

A —^ p. (5.3) 


In this limit, the Salpeter equation (|3.13|) requires that 

E = Ao. 

Using (|5.2|), we hnd 

= (ar')^ + ~ (ar')^ + (^)^. 


(5.4) 


(5.5) 


The leading trajectory is given by the lowest energy for a hxed angular momentum ^ =0, 


so 

L = ar . 

With the circular orbit condition, p = L/r', this gives 

= 2aL, 

and the Regge slope becomes 

£2 2a ■ 


(5.6) 


(5.7) 


(5.8) 


We note that this is also the result given by the Dirac equation with scalar conhnement 
T0| in a heavy-light meson and that time component vector conhnement gives a Regge slope 


of 


OL 


vector 


4a' 


(5.9) 


As demonstrated analytically in |]^ , scalar conhnement with the Salpeter equation does 
not give the desired linear trajectories. 
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VI. CONCLUSION 


Our starting point here has been Buchmiiller’s physical picture relating the static chro- 
moelectric flux tube and the resulting pure Thomas spin-orbit conflnement interaction in a 
meson, though in this paper we make use only of the result that the chromomagnetic held 
vanishes at the quark position. When this assumption is expressed in the laboratory frame, 
its full content becomes more apparent. Once the lab frame four vector potential is found 
by Lorentz transforming a purely electric potential in the instantaneous quark rest frame to 
the laboratory frame, a complete relativistic effective Hamiltonian can be directly computed. 
Not only is the expected Thomas spin-orbit term found, but all other relativistic corrections 
can be evaluated. 

The principal virtue of the model presented here is that it has very few assumptions 
and further that these assumptions are strongly motivated by QCD. We assume that the 
only held acting on the quark is a constant radial electric held in its rest frame and that 
the magnetic held at the quark vanishes in its rest frame. From these assumptions follow 
the Thomas type spin-orbit interaction that is strongly favored by lattice simulations and 
phenomenological considerations, as well as linear Regge trajectories. An additional virtue 
of the model is its simplicity. It is essentially a potential model reducing to the solution of 
a radial differential equation in the general case. 

For reason of simplicity and clarity of exposition, we have considered one quark to be 
fixed. Just as for all potential models, our results carry over directly to the two body case 
with suitable introduction of center of momentum and relative coordinates. For a realistic 
meson we must also include a short range vector (Coulomb) interaction. In the heavy-light 
limit this short range term will result in a spin-orbit interaction opposite in sign to the 
Thomas term. 
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